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ON THE STABILITY AND POINTING OF AN ATTACHED 
DOUBLE -GIMBAL EXPERIMENT PACKAGE 

By Nelson J. Groom, John D. Shaughnessy, 
and Vilas D. Nene* 

Langley Research Center 

SUMMARY 

The pointing capability of a double-gimbal experiment isolation and control system 
of the Apollo Telescope Mount (ATM) type is investigated. Three composite structural 
models of an experiment package connected through frictionless gimbals to a carrier 
vehicle (including a rigid-package — rigid-carrier model, a rigid-package — flexible - 
carrier model, and a flexible -package — flexible -carrier model) are used for this inves- 
tigation, and a linear model of the nominal ATM gimbal control system is considered. 
Contributions of the primary carrier control system to experiment pointing are neglected. 
A linear stability analysis is performed to verify stability of the experiment control sys- 
tem with nominal gains. Transform techniques are then used to compute pointing errors 
onboard the experiment package due to random crew motions input into the carrier vehi- 
cle. Results of the investigation indicate that there is no stability problem due to flex- 
ible coupling of the experiment package and the carrier, and that the ATM-type double- 
gimbal experiment isolation and control system is capable of pointing accuracies better 
than 0.1 arc second. 


INTRODUCTION 

By placing telescopes in orbit above the atmosphere, several severe limitations 
that have hampered ground-based astronomy are removed (opacity of the atmosphere in 
certain wavelength regions and scattering of light from scattering centers), and by having 
the telescope station manned, maximum information can be obtained (ref. 1). In order 
to achieve the necessary pointing accuracy for astronomy experiments on a manned vehi- 
cle, however, some method of isolating the experiments from crew -motion-induced dis- 
turbances is usually necessary. One method of isolation and control is to attach the 
experiment package to the main vehicle through two gimbals with the center of mass of 
the package being placed as close as possible to the intersection of the gimbal axes. The 
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package is then controlled with respect to the main vehicle by gimbal torquers. This 
method will be implemented with the Apollo Telescope Mount (ATM), which will be one 
of the experiments on the Skylab mission (ref. 2). The ATM houses solar astronomy 
experiments and has a specified pointing accuracy of 2.5 arc seconds (defined here as 
the accuracy of maintaining a certain acquired attitude as opposed to the accuracy of 
acquiring the attitude). An ATM-type system is also being considered for use with 
advanced stellar telescopes which require a pointing accuracy of 0.1 arc second. 

The ATM double -gimbal concept has two inherent sources of error which limit the 
pointing accuracy that can be achieved. One of these sources results from the fact that 
the center of mass of the experiment package cannot be placed exactly at the intersection 
of the gimbal axes. Consequently, translational acceleration and rotational motion about 
the pointing axis of the package produces an inertia torque which acts as a disturbance. 
The other source is due to the design itself; since there are only two degrees of freedom 
(rotation about the gimbal axes), the system provides essentially no isolation from trans- 
lational inputs. Thus, disturbance in the carrier vehicle could cause the experiment 
package to be subjected to considerable vibration. If the experiment package were per- 
fectly rigid, then the center -of -mass offset would be the only source of error. The 
experiment package, of course, is not perfectly rigid, so bending of the experiment pack- 
age due to translational and rotational motion (about the pointing axis) of the gimbal 
attachment points must be considered as a source of error. 

The ATM -type system has been investigated for a nonplanar flexible carrier and 
rigid experiment package with nonlinear control (ref. 3), and a flexible model of the 
experiment package has been analyzed by using a planar model attached to a rigid sup- 
port with linear control (in an investigation conducted by the Space Support Division of 
Sperry Rand Corporation under contract to George C. Marshall Space Flight Center). 

The present investigation considers the pointing accuracy of a linear model of the 
ATM-type system for three nonplanar composite structural models ((1) a rigid-package — 
rigid-carrier model, (2) a rigid-package — flexible -carrier model, and (3) a flexible- 
package — flexible -carrier model) to determine what effect the basic design limitations of 
the ATM-type system might have on achieving a pointing accuracy of 0. 1 arc second in 
the presence of random crew motions at various points in the carrier vehicle. These 
three models were examined to determine, respectively, (1) the effects of package center - 
of-mass offset, (2) the combined effects of package center-of-mass offset and carrier 
flexibility, and (3) the combined effects of package center-of-mass offset, package flex- 
ibility, and carrier flexibility. 

The method of computing the vibration modes of the coupled carrier — experiment- 
package structure is described in appendix A by William C. Walton, Jr., of the Langley 
Research Center. 


2 



Crew motions are represented by a filtered random number sequence (with Gaussian 
distribution and negligible mean) which describes continuous motion of a single astronaut 
performing different tasks. Pointing errors were computed by using transform tech- 
niques. Errors were computed on a 3-sigma basis at the fine sun sensor and each exper- 
iment location in pitch and yaw. Only pitch and yaw were used since the experiment con- 
trol system exerts no active control in roll. 

SYMBOLS 


[A] coefficient matrix 

{b} = [a]{t| 

[CONN] matrix defining the connection between carrier and telescope 
[CUP] coupling matrix for carrier and telescope modes 
{d} translations and rotations for coupled system 

{D[ Laplace transform of {d[ 

{DCf, [DC] deformation vectors used in calculation of carrier modes 

IDT}, [DTj deformation vectors used in calculation of telescope modes 

{DRAPC} deflection -rotation vector for telescope attachment point in carrier 
coordinates 

]DRAPT[ deflection-rotation vector for telescope attachment point in telescope 
coordinates 

e base of natural system of logarithms, 2.71828 

]e[ vector of total generalized force for coupled system 

[ECC] telescope mass-point coordinate matrix 

{EOT}, [EDTj elastic deformation vectors used in calculation of telescope modes 
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\e} 

f 

f(t) 

f(k) 


FG) 

F(o>) 

F(s) 

[gmc] 

[gmt] 

h(t) 

H 1> H 2 

H s 

H(oj) 

H(s) 

i=fi 

[I] 

Im 


Laplace transform of |e} 
frequency variable 

input to system, a continuous function of time 

an array obtained by sampling function f(t) 

folding frequency 

fundamental frequency 

sampling frequency 

discrete Fourier transform of f(k) 

continuous Fourier transform of f(t) 

Laplace transform of input to crew-motion filters 
diagonal matrix of generalized carrier masses 
diagonal matrix of generalized telescope masses 
impulse response of system, a continuous function of time 
transfer function blocks in system feedback loop 
fine -sun -sensor gain 
continuous Fourier transform of h(t) 
generalized transfer function in s domain 

identity matrix 
imaginary part of 
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j,k,m,n integers 

K represents a given transfer function block in system forward loop (with 

subscripts 1 to 12) 

K a K m torque -motor gain 

K rg rate -gyro gain 

[KC] symmetric stiffness matrix used to compute carrier modes 

[KT] symmetric stiffness matrix used to compute telescope modes 

L Lagrangian 

[m] = [I]s2 + 2(0.01)[w n ]s + [w n 2 ] 

[MC] symmetric mass matrix used to compute carrier modes 

[mc] carrier mass matrix [MC] with telescope mass deleted 

[MT] symmetric mass matrix used to compute telescope modes 

N number of points in sequences f(k), r(m), and so forth 

|qj vector of modal deflections for coupled system 

{Q} Laplace transform of jq} 

{QC}, LQCj vectors of carrier mode coordinates 
{QT}, LQTj vectors of telescope mode coordinates 
r(m) an array obtained by sampling function r(t) 

r(t) output of system, a continuous function of time 

R(m) discrete Fourier transform of r(m) 
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RM 


R(s) 


Re 


continuous Fourier transform of r(t) 


Laplace transform of linear-system output 


real part of 


Laplace operator 


time variable 


{t} 

[TAP] 

[U] 

[UD] 

[UF] 


periodicity of f(t), if f(t) is periodic, or period of interest 


transfer functions obtained by dividing given e by given input 


6X6 telescope -attachment -point transformation matrix 


modal matrix for coupled system 


modal participation matrix for coupled system 


modal distribution matrix for coupled system 


|u<3c] carrier mode shapes (eigenvectors) for coupled system 
]UQt| telescope mode shapes (eigenvectors) for coupled system 
[UQC], [UQT] matrices of eigenvectors 


[UC] 


matrix of carrier modes 


[UT] matrix of telescope modes 

[wSQc] diagonal matrix of squares of carrier-mode natural frequencies 
[wsqtJ diagonal matrix of squares of telescope-mode natural frequencies 
X c ,Y c ,Zc carrier axes 


Yrj,, Zpji 


telescope axes 
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sampling interval in frequency domain 


sampling interval in time domain 

error point in block diagram representation of system 
(with subscripts 1 to 20) 

command input to Yrp-axis fine sun sensor 

command input to Z^-axis fine sun sensor 

pointing error at fine -sun -sensor location about Yrp-axis 

pointing error at fine-sun-sensor location about Zrp-axis 

rate -gyro damping ratio 

damping ratios in crew-motion filter equations 

standard deviation 

torque -motor time constant 

fine-sun-sensor time constant 

time constants in crew-motion filter equations 

frequency variable 

natural frequencies in crew-motion filter equations 

folding frequency 

rate -gyro natural frequency 


natural frequencies for coupled-system vibration modes 



Notation: 


[ ] 

rectangular matrix 

[ f 

transpose of [ ] 

[ r 1 

inverse of [ ] 

[ i ] 

partitioned matrix 

E 3 

diagonal matrix 

i i 

column vector 

i j 

row vector 

HI i 

partitioned vectors 


A single dot over a symbol denotes a first derivative with respect to time. Double 
dots denote a second derivative. 

ABBREVIATIONS 

AM airlock module 

ATM Apollo Telescope Mount 

CFT continuous Fourier transform 

CMG control-moment gyroscope 

CSM command and service module 

c.m. center of mass 

DFT discrete Fourier transform 

FFT fast Fourier transform 
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LM 

lunar module 

MDA 

multiple docking adapter 

OWS 

orbital workshop 


DESCRIPTION OF STRUCTURAL MODELS 

There are three structural models used in this paper. Each represents a mathe- 
matical model of the ATM experiment package connected through frictionless gimbals to 
a preliminary version of the Sky lab space station. 

One model, henceforth called the flexible-body model, consists of a flexible ATM 
package connected to a flexible model of the Skylab carrier vehicle. The mathematical 
model was formulated from independent structural models for the Skylab and ATM spar 
canister by using the method given in appendix A. The independent structural models 
were developed under separate contracts to George C. Marshall Space Flight Center. The 
resulting composite model consisted of 107 flexible-body modes, with frequencies ranging 
from 0.135 Hz to 38.443 Hz, and eight rigid-body modes (six for the Skylab and two for the 
ATM package). 

The second model, henceforth called the rigid-ATM model, consists of a rigid ATM 
package connected to a flexible model of the Skylab carrier vehicle. This model was 
developed by using the data from the flexible -body formulation except that the flexible 
ATM package modes were omitted. The resulting composite model had 93 flexible -body 
modes, with frequencies ranging from 0.135 Hz to 22.777 Hz, and the same eight rigid- 
body modes as in the model just described. 

A third model, henceforth called the rigid-body model, represents a rigid ATM pack- 
age connected to a rigid model of the Skylab carrier vehicle. This mathematical model 
was obtained by dropping all the flexible-body modes. The resulting composite model had 
the same eight rigid-body modes as in the other two models. 

The Skylab configuration considered, shown in figure 1, is composed of a Saturn IV-B 
orbital workshop (OWS), a multiple docking adapter (MDA), a command and service module 
(CSM), and a modified lunar module (LM) with a rack which provides attachment points for 
solar panels, control -moment-gyroscope (CMG) hardware, and the ATM gimbal suspension 
system. 

The example ATM configuration shown in figure 2 is composed of a double -gimbal 
suspension system with frictionless zero stiffness bearings, an outer shell canister, and 
a cruciform spar to which eight solar astronomy experiments, four rate -gyro packages, 
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one two-axis sun sensor, and other equipment are attached. At each gimbal point, a 
direct-current (dc) torquer is used to provide control torques for the ATM package. 

The mathematical representation of the structural models is based on the assump- 
tion usually made (see refs. 3 and 4, for example) that structural damping, which is not 
included in the modal analysis, can be added in an uncoupled sense in modal coordinates. 
A value of 1 percent of critical damping is assumed to act on each flexible mode. The 
use of this value is consistent with previous analyses (for example, refs. 3 and 5). This 
formulation results in the equations of motion given as 

{q}+ 2(0.01) [o> n ]{q} + [w n 2 ]fa} = [UF] T {Ef (1) 

Here, jq} is the 115 -component vector of modal displacements, £o> n J is the diagonal 

matrix of modal frequencies, and [UF]^ is the matrix of modal data that distributes the 
force vector {e} over the 115 modes and is found by transposing the matrix [UF] 
formed from the modal matrix of the coupled carrier -experiment system [U] (see 
eq. (A 13)) by retaining only those rows corresponding to force or torque application loca- 
tions and directions. The resulting motion in carrier coordinates is given by the relation 

|d}=[UD]{qf (2) 

where [UD] is the matrix formed from [U] by retaining those rows that correspond to 
the motion of the mass point of interest. 

DESCRIPTION OF CONTROL SYSTEM 

The experiment control system used in this paper is a linear model of the ATM sys- 
tem with the roll control system fixed (the roll system is manually controlled prior to 
experiment initiation). The system consists of a two-degree-of-freedom vernier gimbal 
system driven by two dc torque motors per axis with error signals derived from a two- 
axis fine sun sensor and rate gyros mounted on the experiment package. A block diagram 
of the system used is presented in figure 3, and table I gives the parameters used. 

Since the CMG system, which provides Skylab spacecraft stabilization, is designed 
to control rigid-body motion resulting from low-frequency, low -magnitude torques instead 
of crew-motion disturbances (ref. 3), no primary CMG control system was included. 
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CREW-MOTION MODEL 


Previously, deterministic approximations of typical astronaut movements were 
employed to represent crew -motion disturbances (see appendix II of ref. 6). Recent 
studies have developed linear filters that can be used in conjunction with a unity-power- 
spectral-density noise source with Gaussian distribution and negligible mean to represent 
complete activities such as the performance of a given experiment. These activities 
combine most of the typical astronaut movements described in reference 6. 

The crew -motion input disturbances for this paper were obtained by passing a 
computer -generated pseudo-random number sequence with a unity power spectral den- 
sity and Gaussian distribution with negligible mean through appropriate linear filters. 
Three separate activities were chosen as being representative. These activities con- 
sisted of two manual console operations and one compressive walking activity. All three 
are modeled as single -astronaut motion. The resultant filters have a general transfer 
function H(s) given by 


H(s) = 


t 1 s ( t 2 s+ X ) 


f 2 
s z + 


2 P 






S 2 + 


2 P 2 S 2 s + 


5 o 2 


(3) 


The values of the parameters in equation (3) are given in table II for each crew activity. 
The two manual console operations are restrained activities. That is, the astronaut is 
restrained (by harnesses or other restraining devices) from large translational move- 
ments of his center of mass. These two activities (Run 1 and Run 4) represent an astro- 
naut performing the S 056 and S 082 A ATM experiments (these experiments are defined 
in ref. 1). The compressive walking activity is an unrestrained activity which consists 
of lateral translation between structures where the astronaut can support himself with 
hands and feet simultaneously. The filters and parameters were developed by Martin 
Marietta Corporation under contract to George C. Marshall Space Flight Center. The 
method used to develop the filter transfer functions is given in reference 7. The crew 
motions were input into two different locations in the carrier vehicle: the AM/MDA 
intersection and the OWS crew location (fig. 1). All three activities were investigated 
for the first location, whereas only compressive walking was assumed for the OWS crew 
location. 


METHOD OF ANALYSIS 

The experiment -package pointing errors were computed by convolving the transfer 
function between the input to the crew-motion filters and the fine -sun -sensor outputs 
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(F(s)/e 6 and F(s)/ej 2 i- n fig- sj with a random number sequence. The random number 
sequence had a unity power spectral density and a Gaussian amplitude distribution with 
negligible mean. The frequency -domain transfer -function equations are developed in 
appendix B. The convolution was performed by using discrete Fourier transform (DFT) 
techniques and consisted of (1) computing the numerical value of the transfer function of 
interest at discrete values of jo>, (2) taking the forward DFT of the random number 
sequence representing the input, (3) multiplying point by point the sequence representing 
the numerical value of the transfer function and the DFT of the input, and (4) taking the 
inverse DFT of the product to obtain a sequence representing the error time history of 
interest. The analytical details of this method are given in appendix C. Table HI gives 
the parameters used in the actual computation of error time histories. The 3a values of 
pointing errors were computed from these time histories. All forward and inverse 
transforms were computed by the fast Fourier transform (FFT) algorithm which is some- 
times called the Cooley-Tukey algorithm (ref. 8). All computations were performed on 
a Control Data 6600 computer system. 

RESULTS AND DISCUSSION 
Stability 

A linear stability analysis was performed to verify that the control system being 
considered was stable with the nominal ATM torque -motor gain (2K a K m ). The stability 
analysis was conducted by applying the Nyquist stability criterion to the open-loop trans- 
fer function of each loop ^ C y/ e 6 0 C z/ e l2 in wi -th °ther loop closed. 

The upper limit of 2K a K m for stability was found to be 136.9 N-m/volt. It was expected 
that including the flexible model of the carrier in the stability analysis instead of assuming 
a rigid immovable base would result in a stable gain much lower than the nominal, espe- 
cially since the nominal ATM control loop contains no compensation. However, as the 
results show, stability is not a problem for the model analyzed and the nominal gain 
could, in fact, be increased. 


Rigid- Body Model 

Results obtained by subjecting the rigid-body model to crew-motion disturbances 
are presented in table IV. Since the total angular errors at different points of the exper- 
iment package are the same for the rigid-body case, only the fine -sun -sensor errors are 
given. These results represent errors due entirely to the center -of-mass offset since 
no flexibility is involved. The center-of-mass offset for the experiment -package model 
(in experiment -package coordinates) was computed from the ATM mass data as 
(-0.000914, 0.02275, -0.02159) in meters. This amount of offset can be considered to be 
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excessive since it is possible to reduce the offset to 0.00127 meter along each axis by 
careful balancing (ref. 3). 


Rigid-ATM Model 

Results obtained by subjecting the rigid-ATM model to crew-motion disturbances 
are presented in table V. Again, only the fine-sun-sensor errors are given because the 
experiment package is rigid. These results represent combined effects of experiment- 
package center-of-mass offset and carrier -vehicle flexibility. In some cases, the addi- 
tion of the carrier flexibility caused an increase in the 3a pointing error of almost an 
order of magnitude. The maximum error, however, was still well below 0.1 arc second. 

Flexible -Body Model 

Results obtained by subjecting the flexible-body model to crew-motion disturbances 
are presented in table VI. Again, only the fine -sun -sensor errors are given because it 
was found that, although the total angular errors at different points of the experiment 
package were not exactly equal, the differences were insignificant. These results repre- 
sent errors due to the combined effects of experiment -package center-of-mass offset, 
experiment-package flexibility, and carrier -vehicle flexibility. In some cases, the 
experiment -package flexibility caused an increase in the 3a pointing error by a factor of 
approximately 3 over the rigid-ATM-model errors. Thus, not only bending of the carrier 
but also bending of the experiment package appears to be a significant contributor to 
pointing error. Also, since the maximum 3a error was approximately 0.02 arc second, 
the basic design limitations of the ATM-type system appear to present no obstacle in 
achieving a pointing accuracy of 0.1 arc second. 

Multiple Crew Motion 

The errors for each structural model are computed for single -astronaut motion. 
Since there will be more than one crew member, it is likely that there will be more than 
one astronaut disturbance at any given time. Since the systems analyzed are linear 
systems, a worst-case estimate of errors due to multiple crew motion can be made by 
simply adding 3a errors for different single -astronaut disturbances. 

Gain Sensitivity 

In order to determine the sensitivity of the gimbal control loops to changes in gain, 
runs with different gains were made by using one carrier location and one crew -motion 
input. The carrier location was the AM/MDA intersection and the crew motion was 
Console Operation Run 1. A plot of 3a error as a function of gain is shown in figure 4 
for the fine -sun- sensor Y^-axis and the fine -sun -sensor Zrp-axis. The sensitivity of the 
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pointing error of the basic ATM -type gimbal control loop to changes in torque-motor 
gain suggests that the application of optimal control concepts could result in a significant 
reduction of overall pointing error. 

CONCLUDING REMARKS 

In order to examine the inherent design limitations of a double -gimbal experiment 
isolation and control system of the Apollo Telescope Mount (ATM) type, three mathemat- 
ical structural models of an experiment package (ATM) connected through frictionless 
gimbals to a flexible carrier vehicle (Skylab) were utilized. The three models consisted 
of the following: 

(1) The rigid-body model of the experiment package connected to the rigid-body 

model of the carrier vehicle 

(2) The rigid-body model of the experiment package connected to the flexible model 

of the carrier vehicle 

(3) The flexible model of the experiment package connected to the flexible model of 

the carrier vehicle 

A linear model of the nominal ATM gimbal control system with no sensor, actuator, 
or bearing nonlinearities was used, and no control-moment-gyroscope control system for 
the carrier vehicle was included. The stability of the ATM gimbal control system was 
investigated on the model of the flexible experiment package connected to the flexible 
carrier vehicle since, on this model, the actuators would be torquing against a flexible 
structure instead of a fixed rigid base as is normally assumed. The models were sub- 
jected to continuous single -astronaut crew motion, represented by a filtered random num- 
ber sequence with Gaussian distribution and negligible mean, and pointing errors were 
computed on a 3-sigma basis about the two gimbal axes of the experiment package at the 
fine-sun-sensor location. 

The results of analyses performed on the composite models indicate the following: 

1. No stability problem is expected to result from flexible coupling of the experiment 
package and carrier vehicle. 

2. There is enough coupling between flexible-experiment-package modes and 
flexible -carrier modes that bending of the experiment package must be considered a sig- 
nificant contributor to experiment-package pointing error. 

3. The ATM double-gimbal experiment isolation and control system design is inher- 
ently capable of providing pointing accuracy better than 0.1 arc second. 
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4. The application of optimal control concepts to the basic ATM -type gimbal control 
loop could result in a significant reduction of overall pointing error. 

In practice, sensor, actuator, and bearing anomalies, which were neglected for this 
study, could be limiting factors as far as achieving 0.1-arc-second pointing accuracy in 
the presence of crew disturbances. A nonlinear study of the double-gimbal system, 
utilizing the composite models developed in this report, would be required to determine 
what effect these anomalies would have. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 24, 1972. 
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APPENDIX A 


METHOD FOR COMPUTING VIBRATION MODES OF 
THE CARRIER-EXPERIMENT SYSTEM 

By William C. Walton, Jr. 

Langley Research Center 

This appendix explains the method used in the present paper to obtain a mathemat- 
ical model of a manned carrier vehicle (Sky lab) carrying a large telescope (ATM). The 
math model is a modal model consisting of computed natural vibration modes of struc- 
tural components. The method has been specially devised to use vibration modes pre- 
viously computed for the carrier vehicle with the ATM experiment package considered 
as a rigid attached mass and vibration modes for the ATM considered as an elastic body 
attached to a rigid base. Special features of the application of the method to the coupled 
carrier — experiment -package structure are given. 

The following assumptions are made: 

(1) Computed vibration modes and natural frequencies are available for the carrier 
structure with the experiment package included as a rigid body attached to the carrier 
through locked gimbals. These modes will henceforth be referred to as carrier modes. 

(2) Computed vibration modes and natural frequencies are available for the experi- 
ment package considered as an elastic body gimbaled to a fixed base. These modes will 
henceforth be referred to as ATM modes. 

(3) Lumped-mass analyses were used in obtaining both carrier and ATM modes, 
and the computed mode shapes are expressed as deflections and rotations at the mass 
points (see ref. 4, ch. 5, for treatment of vibrations of lumped-mass systems). 

(4) The masses and mass-point locations utilized in the computation of the ATM 
modes are available. 

(5) Information necessary to define the connection between the carrier and ATM 
package is available. 

The Lagrangian for the coupled system may now be developed. The vector {DC} 
is defined as the deformation vector used in the calculation of the carrier modes. It is 
assumed that 

{DC} = [UC]{QC{ (Al) 
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APPENDIX A - Continued 

where the columns of [UC] are carrier modes and {QC} is the vector of carrier mode 
coordinates. The matrix [UC] is generally rectangular. The vector IDT} is defined 
as the deformation vector used in the calculation of ATM modes. It is assumed that 

|dt[ = [CONN]|DC} + [UT]]QT} (A2) 

where the columns of [UT] are ATM modes, |QT[ is a vector of ATM mode coordi- 
nates, and [CONN] is a matrix to be developed subsequently such that the term 
[CONN]|DC[ represents the motion of the ATM package as if it were rigidly attached to 
the carrier. 

The modal analyses which produce carrier modes and ATM modes are based on 
structural idealizations which are presumed to imply stiffness and mass matrices. 
Although some of these matrices are not needed in computation, their use is convenient 
in derivations and they are defined here for that purpose. 

Let 


{EOT} = [UT]|QT| (A3) 

Then, the Lagrangian for the coupled system may be written as 

l = |LdcJ[kc]|dc[ + |ledt][kt]{edt} - |[dcJ[mc]|dc} - |[dtJ[mt]|dt}- (A4) 

Matrices [KC] and [MC] are, respectively, the symmetric stiffness and mass matrices 
used to compute carrier modes; [KT] and [MT] are, respectively, the symmetric 
stiffness and mass matrices used to compute ATM modes; and £mc] is the carrier mass 
matrix [MC] with the ATM mass deleted. 

By referring to equation (A2), equation (A4) may be rewritten as 

l = I(-[dcJ[mc]{dc} - [pic] [ c onn] t [mt] [c onn]|d'c | 

- |d'cJ [c onn] t [mt]|edt| - |edtJ [mt] [conn] |dc } 

- [edtJ[mt]|edt| + Ldc][kc]{dc} + Ledt][kt]{edt^ (A5) 

The second term within the parentheses on the right-hand side of equation (A5) represents 
the kinetic energy of a motion of the ATM as if the ATM were rigidly attached to the car- 
rier. Therefore, by considering the previously given definition of J~McJ, the first two 
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(A6) 


(A7) 


terms within the parentheses may be combined to give 

l = |^dcJ[kc]{dc}+ |edt|[kt]{edt} 

- |pc] [mc]|dc| - [dcJ [conn] t [mt]{edt| 

- [edtJ[mt][conn]|d’c| - [edtJ[mt]|edt}J 

By using equations (Al) and (A3), equation (A6) can be written as 

L = |(lQCj[UC] T [KC][UC]{QC} + [qtJ [ut] t [kt] [ut] {qt } 

- |qc] [uc] t [mc] [uc]{qc| - [qcJ[uc] t [conn] t [mt] [ut]{qt| 

- [qtJ[ut] t [mt] t [conn][uc]{qc} - [qt] [UT] T [MT][UT]{qt} j 

It is now noted that certain of the expressions which occur in equation (A7) may be 
written as follows: 

[UC] T [MC][UC] = [gMC^] (A8a) 

[UT] T [MT][UT] = [gMt] (A 8b) 

[UC] T [KC][UC] = ^GMcJ^WSQC^ (A8c) 

[UT] T [KT][UT] = [GMTj [wSQT^| (A8d) 

where ^GMcJ, ^ GMtJ, |^WSQC^, and ^WSQT^] are diagonal. The diagonal ele- 
ments of J^GMC^j and [gMtJ are the generalized masses associated with the carrier 
modes and the ATM modes, respectively. The diagonal elements of ^WSQC^ and 
[wsqtJ are the squares of the natural frequencies associated with the respective 
modes. The coupling matrix [CUP] is defined by the equation 
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[CUP] = [uc] t [conn][mt][ut] (A9) 

By using equations (A3), (A8), and (A9) in equation (A7), the Lagrangian becomes 


= |[qc | qtJ 


[gmc] [wsqc] 

0 

1 

o 

[gMt] [\VSQt] 


QC 

QT 


QC1QT 


[gmcJ ; 
1 

[CUP] 

[cup] t i 

i 

1 

1 — 71 
H 
2 
o 

\L 1 


QC 

QT 


(A10) 


Equation (A10) leads to the eigenvalue -eigenvector problem 


[gmc] [wsqc ^ 

]j 0 

- 0) 2 

[gmc] 

0 

i [gmt] [wsqt] 


[CUP] T 



UQC | 

1 = 0 (All) 


UQT 


The eigenvalues w n 2 of equation (All) determine the natural frequencies of the coupled 

carrier -experiment system and the associated eigenvectors [uQCjUQt] are the mode 
shapes in terms of modal coordinates. To express the mode shapes in terms of physical 
displacements and rotations, the following equation derived from equations (Al) and (A2) 
may be used: 


DC 


1 

[UC] i o 


UQC 

DT 


i 

[CONN][UC] i [UT] 

1 

1 _ 


UQT 


(A12) 


The model matrix [u] for the coupled carrier -experiment system which is referred to 
in the text is thus calculated from the equation 
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M- 


[uc] 

„ - J 

0 

UQC 

[conn] [uc] 

[UT] 

UQT 


(A13) 


where the columns of the matrix on the far right constitute a complete set of eigenvectors 
of equation (All). 

The matrix [CONN] is derived on the assumption that the ATM package and the 
carrier vehicle are connected at only one point and that the common point is a mass point 
of the carrier and also a mass point of the ATM experiment package. Thus, the equation 
of constraint for locked gimbals can be expressed as 

|DRAPT| = [TAP] |DRAPC[ (A 14) 

where the elements of the vectors |DRAPT| and |DRAPC] are deflections and rota- 
tions of the attachment point referred to the ATM coordinate system and to the carrier 
coordinate system, respectively, and where [TAP] is a 6 x 6 transformation matrix 
depending only on the relative orientations of the two coordinate systems. Since the 
attachment point coincides with mass points of both the ATM and the carrier, the elements 
of {dRAPT] are six consecutive elements of the vector |DT}, and the elements of 
■{DRAPC] are six consecutive elements of |DC}. If it is assumed that the ATM is rigid, 
then the vector |DT} is related to the attachment -point displacements and rotations by 
an equation of the form 

|DT[ = [ECC] {DRAPT} 

= [ECC][TAP]{DRAPC[ (A 15) 

The matrix [ECC] depends only on the coordinates of the mass points of the ATM 
referred to the ATM coordinate system. Thus, the matrix [CONN] has the form 

CONN = [O ! [ECC] [TAP] ! oj (A 16) 

where the finite columns correspond to the six elements in {DC] which are identical to 
the elements of {DRAPC}. 
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FORMULATION OF SYSTEM-TRANSFER- FUNCTION EQUATIONS 

Modal Equations 

The basic uncoupled modal differential equations for the composite model are given 
by equation (1) as 


{q}+ 2(0.01) [w n ]{q}+ [w n 2 ]{q}= [UF] T {E| 


If the Laplace transform of this equation is taken, it becomes 


([I]s 2 + 2(0.01)[co n ]s + [w n 2 ] J{Q} = [UF] T {E} 


(Bl) 


where [I] is the identity matrix, s is the Laplace operator, {Q} is the Laplace 
transform of jq}, and {e} is the Laplace transform of |e|. For simplification, the 


diagonal matrix 



is defined as 


*M]= [I]s 2 + 2(0.01)[ w n Js + [ w n 2 ] 


(B2) 


Then, 

{Q}=[mJ[UF] T {e} (B3) 

If |Df represents the Laplace transform of the modal displacements in carrier coordi- 
nates found by transforming equation (2), then 


|D}= [UD]{Qf 

Finally, substituting from equation (B3) results in 


(B4) 


|d}= [ud][m J’ 1 [uf] t {e} 


(B5) 



APPENDIX B - Continued 


System Transfer Function 

A block diagram of the experiment -package Y T - and Z T ~axis servo loops and one 
general set of crew-motion filters is given in figure 3. The crew-motion filters are the 
blocks labeled Ky through Kj 2 - First, the transfer -function equations for the crew- 
motion input, labeled F(s), will be developed and then the transfer-function equations for 
the stability analysis will be developed from these equations. In figure 3, eg, eg, e 
and e 2 Q correspond to {D} in equation (B5) and eg, e^, eg, e jq, and e ^ through 
e jg correspond to {1} in the same equation. With reference to figure 3, the matrix 
equation that includes each e (with 6> cy and 0 CZ set to zero and with input F(s)) is 


F(s) 
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APPENDIX B - Continued 

Note that in figure 3 each torque motor has two outputs. These outputs represent a torque 
input into the experiment package and a reaction torque input into the carrier. Since the 
experiment -package roll ring is considered attached to the carrier in a rigid plane (see 
appendix A), the reaction torques for one axis go into one location in the carrier. Dividing 
through the matrix equation (B6) by F(s), and then rearranging, results in 


l 
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> (B7) 


where a given element in {t} represents a transfer function with F(s) as the input 

( for example, T 9 = V Equation (B7) has the form 
\ Z F(s)/ 

= {B} 
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The solution for {t} for a given value of s then is 

= <B8) 

and T 6 , T 12 , T 5 , and Til represent, respectively, the transfer function between the 
crew-motion input and the angle about the Yrp-axis at the fine-sun-sensor location, the 
angle about the Z T ~axis at the fine -sun -sensor location, the rate about the Y^-axis at the 
rate-gyro location, and the rate about the Z^-axis at the rate-gyro location. The equa- 
tions for these four functions are used in finding the time response. The method for 
finding the response is presented in appendix C. 

By changing the crew-motion input equations slightly, the experiment -package 
Yrp-axis sun-sensor open-loop transfer function and Zrp-axis sun-sensor open-loop trans- 
fer function can be found. To find the Y T -axis sun-sensor open-loop transfer function, 
set 0 C y and F(s) to zero and change the equation for e 2 as follows: 

e 2 = 0 cy K l " e 5 H l 

Dividing by 0 C y and transposing gives 
*2 + * 5 "! - K 1 

With this change, equation (B7) can be suitably modified. The Zrp-axis sun-sensor open- 
loop transfer function can be similarly obtained. 
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APPENDIX C 


OBTAINING LINEAR-SYSTEM RESPONSE BY USING 
THE DISCRETE FOURIER TRANSFORM 

A linear system can be represented as 


f(t) 


Random signal 
input 


H(s) = R(s)/F(s) 

System transfer 
function 


r(t) 


Pointing error 
output 


The input f (t) to the system and the transfer function H(s) of the system are assumed 
to be known. This appendix explains how to obtain the output r(t) of the system by 
using discrete Fourier transform techniques. 


Continuous and Discrete Fourier Transforms 

If f(t) is defined for a period 0 st ^T and is zero elsewhere, the continuous 
Fourier transform F(w) can be obtained as 

F(co)=f f(t)e~ iwt dt (Cl) 

— oO 

That is, 

T 

F(co) = C f(t)e" iwt dt (C2) 

J 0 

The transform F(w) can be obtained for any frequency from to 00 and then the 

time function f(t) can be obtained as 

f(t) = i-C F(co)e iu,t du> (C3) 

" TT *J _oo 

Now, if F(u>) is obtained only for the frequency range and if Wj is suf- 

ficiently high to make F(w) practically zero for all |o>| > | w f|> then it is possible to 
write 
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f(t) = A- f f F(co)e ia ’ t dw (C4) 

2ir 

The two functions f(t) and F(a>) can be considered as a CFT pair and may be 
denoted by the symbol f(t) — > • F(o>). 

The system output r(t) can be obtained by convolving the impulse response h(t) 
with the input f(t). Therefore, if r(t) R(a>) and h(t) — H(oj), then 

R(o>) = H(w)F(w) (C5) 


and 


r(t) = — f H(w)F(co)e icot da) (C6) 

2tt J-oo 

If H(u>) and F(w) are available for I w ^ then 
r(t) = — H(w)F(w)e iwt da> 

2 TT 

It can be seen that H(cd) can be obtained from the transfer function 
system by replacing s by ico; that is, 

H(w) = H(iw) (C8) 

Thus, if the input f(t) and the transfer function H(s) are known, it is possible to obtain 
the output r(t). 

However, if the system is to be analyzed on a digital computer, one must use the 
finite, discrete version of the Fourier transform (that is, the discrete Fourier transform 
(DFT)), and all the functions or the waveforms have to be sampled and cannot be continuous. 

If the input f(t) is defined for 0 i t g T and sampled at a uniform interval of 
time At, an array f(k), with k = 1,2,. . .,N, will be obtained defining the function f(t). 
Now, analogous to the CFT pair, a DFT pair f(k) F(j) can be defined as 

N -27ri(k-l)(j-l) 

F(j) = N X f(k)e N (j = 1,2, ‘ • -’ N) (C9) 

k=l 


(C7) 

H(s) of the 
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and 


N -2ai(k-l)fl-l) 

ffr)= £ F(j)e N (k = 1,2,. . .,N) (CIO) 

3=1 

where f(k) and F(j) are, in general, complex sequences defining functions in the time 
domain and frequency domain, respectively. 

It can be seen that time and frequency are not explicit in equations (C9) and (CIO). 
However, the index k refers to a sample period number and the index j represents a 
harmonic number. The true time and frequency corresponding to these indices can be 
obtained as 


t k =(k-l)At (Cll) 

and 

fj = G - l)Af (C12) 

where 

Af=( 0 =i (C13) 

At this stage, it should also be noted that the sequences f(k) and F(j) are peri- 
odic over N because of the exponential term in equations (C9) and (CIO). As explained 
subsequently, this property of the DFT has important implications when two time 
sequences are convolved. 

Now, if f(t) is a real time function, the sequence f(k) is real and the sequence 
F(j) has the following properties (see fig. 5): 

(1) The real part of F(j) is symmetric about fj = f 0 , called the "folding 

Ci 

frequency"; that is, 

Re F(j) = Re F(N + 2 - j) (j = 2,3,. . (C14) 

(2) The imaginary part of F(j) is antisymmetric about the folding frequency f^; 
that is, 

Im FG) = -Im F(N + 2 - j) (j = 2,3,. . (C15) 
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The two numbers F(l) and F^ + 1^ in the array F(j) correspond to zero and 
folding frequency, respectively. 

Now, since the sequence F(j) is periodic, one can say that the real part of F(j) 
is an even function and that the imaginary part of F(j) is an odd function. In other 

words, the Fourier coefficients F^ + 2j,. . .,F(N) can be viewed as the negative fre- 
quency harmonics for frequencies between -^fj - f Q ) and -f Q . 

Periodic and Aperiodic Convolutions Using the DFT 

When the functions in the time domain are sampled, the convolution of two sequences 
f ^(k) and f 2 (k), which in the frequency domain is represented by the product of their 

DFT's, takes the following form: 

N N 

FjOFgCj) — ^ f 1 (n)f 2 (j - n + 1) = 1 ^ fjQ - n + l)f 2 (n) (j = 1,2,. . .,N) (C16) 

n=l n=l 

As pointed out earlier, the sequences fj(k), f 2 (k), F^(j), and FgCj) are periodic 

over N; that is, if f(k) — ~ F(j), then F(j) is a periodic function over a frequency band 

0 to f s and describes the DFT of a time function f(k) which is periodic over a period 

0 to T. As a result, the convolution defined by equation (C16) is also periodic. 

Now, if the input f(t) which is sampled over 0 to T is periodic over T, then 
the steady-state output r(t) for a linear system will also be periodic over T and the 
convolution defined in equation (C16) will correctly describe the recurrent transients 
of r(t). However, the use of equation (C16) when f(t) is not periodic over T will be 
incorrect and the sequences have to be suitably modified. 

To illustrate the unwanted periodicity in convolution, assume two sequences f^(k) 
and f 2 (k) as shown in figure 6. The N-point convolution fg(k) of these two sequences 
can be written as 

N N 

f 3 (k) = n Z f l (n)f 2< k " n + ^ = n Z f l (k - n + 1)f 2 (n) (k = 1 ’ 2 ’* * *’ N) (C17) 

n=l n=l 

This computation corresponds to reversing one sequence and sliding it past the 
other in increments as the summing operation of equation (C17) is performed. Let f 2 (k) 
be fixed. Figure 7 shows the sequence f^(k) reversed in time and positioned for com- 
putation of the first point fg(l). The first two points can be written as 
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f 3 (1) = + f l^ f 2^ 2) + f i^- 1 ) f 2 < 3 > 

+ f 1 (-2)f 2 (4) + f 1 (-3)f 2 (5) + f 1 (-4)f 2 (6)| (CIS) 

and 

f 3 (2) = ^i(2)f 2 (l) + f 1 (l)f 2 (2) + f i(0)f 2 (3) 

+ f 1 (-l)f 2 (4) + f 1 (-2)f 2 (5) + f 1 (-3)f 2 (6)| (C19) 

Notice that the negative indices can be interpreted by the following equation because 
of periodicity of the sequences 

f x (-j) = f x CN - j) (j = 0,1,2,. . .,(N-2)) (C20) 

As each point of the reversed sequence f j(k) is shifted out to the right, the corre- 
sponding point is shifted in on the left; hence, unless the sequence being shifted is peri- 
odic over N, the convolution performed by equation (C17) will be erroneous. This situ- 
ation can be easily remedied and an aperiodic convolution can be obtained as follows. If 
two sequences are described by N points, N more points having zero value can be 
added to the sequence to make a total of 2N elements in the sequence. The DFT, the 
inverse DFT, and hence the convolution will then be periodic over 2N, and an N-point 
convolution of the two sequences defined in this manner will give an aperiodic convolu- 
tion of two N-point sequences. This method is illustrated in figures 8 and 9 and described 
in detail in reference 9. 


Computation of the Sampled Output 

A step-by-step procedure has been described below to obtain the system output r(t) 
when the input f(t) and the transfer function H(s) are known: 

(1) Choose a value for N such that the folding frequency fj = — is high enough 

2T 

to describe the transfer function adequately. 

(2) Sample the input f(t) over T at a uniform interval of time ^At = and 

obtain the sequence f(k), with k = 1,2,. . .,N. Extend the sequence by adding zero terms 
to a total of 2N terms. Thus, an array f(k), with k= 1,2,. . .,2N, is obtained such 
that 
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f(k) = fgk - l)At] (k= 1,2,. . .,N) 

f(k) = 0 (k = (N+l),(N+2),. . .,2N) 


(3) Operate the DFT of f(k) and obtain the sequence F(j), with j = 
The increment in frequency corresponding to the two consecutive values of 


then 


*o = 


2T 


1,2,. . .,2N. 

FG) is 


(4) Define the sequence H(m) as follows: 


H(m) = H|i27r(m - l)fjj (m = 1,2,. . .,(N+1)) 

H(N + 1 + j) = conj H(N + 1 - j) G = 1,2, . . .,(N-1)) 


(5) Obtain the array 

R(m) = H(m)F(m) (m = 1,2,. . .,2N) 

(6) Operate the inverse DFT on R(m). 

(7) Form the array r(m) such that 

r(m) = Re R(m) (m = 1,2,. . .,N) 

The array r(m) defines the sampled output r(t) over the period 0 ^ t § T. 
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TABLE I.- NOMINAL PARAMETERS OF THE ATM CONTROL SYSTEM 


Fine-sun-sensor gain, H s , volts/rad 33 000 

Rate -gyro gain, K r g, volt s/rad/ sec 2580 

Torque-motor gain, KgK^, N-m/volt 37.95 

Fine-sun-sensor time constant, r s , sec 0.016 

Torque -motor time constant, r a , sec 0.0032 

Rate -gyro natural frequency, u>g, rad/sec 154.88 

Rate-gyro damping ratio, Pg 0.8 


32 



TABLE II.- CREW-MOTION FILTER PARAMETERS 


Crew motion 

Input 

Input 

axis 

Pi 

W 1 

p 2 

(a) 

*2 

(a) 

T 1 

t 2 



X C 

0.5056 

4.0291 

0 

0 

4.1091 

0 

Console 

operation 

Force 

Y C 

.5565 

2.8065 

0 

0 

3.8489 

0 


z c 

.5986 

4.7634 

0.2357 

31.2563 

2563.12 

0 

Run 1 


x c 

0.3196 

4.4158 

0 

0 

7.27 

0 


Moment 

Y c 

.4683 

4.0384 

0 

0 

12.7861 

0 



z c 

.5589 

4.4132 

0 

0 

2.9666 

0 



x c 

0.6176 

4.8207 

0 ! 

0 

6.2389 

0 


Force 

Y c 

.4660 

3.9640 

0 

0 

5.6195 

0 

Console 

operation 

Run 4 


z c 

.7719 

5.3108 

0 

0 

5.6142 

0 


x c 

0.4846 

4.5882 

0 

0 

13.8367 

0 


Moment 

Y c 

.6750 

4.3448 

0 

0 

20.7737 

0 



z c 

.7691 

5.3496 

0 

0 

5.9213 

0 



x c 

0.2097 

4.7567 

0 

0 

1.2428 

0 


Force 

Y c 

.2601 

14.0260 

0.3831 

2.9387 

375.92 

0.07 

Compressive 

walking 


z c 

.4509 

8.0629 

0 

0 

7.2673 

0 


x c 

0.2509 

1.2705 

0.1231 

13.3204 

931.01 

0 


Moment 

Y c 

.3839 

4.0364 

.3196 

16.9327 

2106.51 

0 



z c 

.3386 

2.7191 

.1850 

13.0889 

476.43 

0 


a If pg and cog are zero, then the term [s? + 2p2&2 s + ^g 2 ) is deleted from 
equation (3). 
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TABLE m.- PARAMETERS FOR COMPUTATION OF 
ERROR TIME HISTORIES 


Number of points in sequences f(k), r(m), etc., N 2048 

Folding frequency, fj, rad/sec 300 

Sampling interval in time domain, At, sec 0.01047 

Period of f(t), T, sec 21.45 


TABLE IV.- RIGID-BODY 3a POINTING ERRORS FOR 
FINE -SUN-SENSOR LOCATION 


Crew motion 

Carrier location 

3a errors, arc sec 

0 ey 

^ez 

Console Operation Run 1 

AM/MDA intersection 

0.3935 x 10“ 3 

0.4160 x 10-3 

Console Operation Run 4 

AM/MDA intersection 

.6273 

.6560 

Compressive walking 

AM/MDA intersection 

.9888 

1.0423 

Compressive walking 

OWS crew location 

.5289 

.5824 
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TABLE V.- RIGID-ATM 3a POINTING ERRORS FOR 
FINE -SUN -SENSOR LOCATION 


Crew motion 

Carrier location 

3a errors, arc sec 

o 

^ez 

Console Operation Run 1 

AM/MDA intersection 

3.8494 x 10-3 

4.2426 X 10“3 

Console Operation Run 4 

AM/MDA intersection 

5.3867 

5.8878 

Compressive walking 

AM/MDA intersection 

7.4736 

7.8396 

Compressive walking 

OWS location 

2.5270 

2.6880 


TABLE VI.- FLEXIBLE -BODY 3a POINTING ERRORS FOR 
FINE -SUN-SENSOR LOCATION 


Crew motion 


3a errors 

, arc sec 

Carrier location 

0 ey 

^ez 


Console Operation Run 1 

AM/MDA intersection 

9.766 X 10-3 

12.672 X 10-3 

Console Operation Run 4 

AM/MDA intersection 

11.971 

15.860 

Compressive walking 

AM/MDA intersection 

15.500 

20.130 

Compressive walking 

OWS crew location 

2.996 

3.218 


* 
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(a) Spar and ring assembly. (Center -of-mass location is given in meters.) 


Figure 2.- ATM configuration. 
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Z T (pitch) 

(b) Gimbal system. 
Figure 2.- Concluded. 















3a pointing errors, arc- sec 



Gain, N-m/volt 

Figure 4.- Variation of total 3a pointing errors with torque-motor gain 

for fine -sun-sensor location. 
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